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ABSTRACT 


Digital imagery acquired by satellite have inherent 
geometrical distortion due to sensor characteristics and to platform 
Variations. In INPE a software system for geometric correction of 
LANDSAT MSS imagery is under development. Such conected imageri will 
be useful for map generation and thematic mapper. Important examples 
are the generation of LANDSAT "Image-charts" for the Amazon region 
and. the possibility of integrating digital satellite imagery into a 
Geographic Information System. 
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1. PROP0SITO do sistema 


1.1 - OBJETIVOS 

Este documento descreve a especificacao e o projeto pre 
liminar de urn sistema para geracao de imagens digital' s de satelite cor 
rigidas geometricamente. Este sistema esta sendo projetado de acordo 
com as mais avancadas tecnicas de engenharia de "software" (Pressman, 
1S82; Yourdon and Constantines, 1979; Page-uones, 1981). Deste modo, 
este documento sera secundado pelo que descrevera o projeto detalhado 
do sistema, o qual incluira diagramas de fluxo e estrutura dos arqui^ 
vos a serem utilizados. Estas tecnicas vem sendo aplicadas em grande 
numero de empresas e institutes na area de Computacao com excelentes 
resultados (Beck and Perkins, 1983). 

0 objttivo do sistema e produzir imagens digitais com 
gidas geometricamente, as quais serao bastante uteis tanto para a ger£ 
cao de cartas quanto para referenciar diretamente a regiao de intere^ 
se a local idades geograficas conhecidas. Exemplos importantes sao a ge 
racao de cartas-imagem LANDSAT para a Regiao Amazonica e a obtencao de 
imagens digitais em projecao desejada para a integracao direta em urn 
Sistema Geografico de Informacao. 

Urn produto adicional a ser gerado pelo sistema serao mo 
saicos digitais, os quais produzirao combinacoes de imagens LANDSAT de 
acordo com as especificacoes desejadas. Urn exemplo de mosaico seria a 
geracao de imagens digitais corrigidas de acordo com a carta UTM pa 
drao na escala 1:250.000. Este documento descreve a especificacao com 
pleta de urn Sistema para Aplicacoes Cartograficas Digitais a ser im 
plantado no DGI, em Cachoeira Paulista, a parti r do suporte computacio 
nal descrito no item seguinte e do suporte log.ico descrito nas Secoes 
seguintes. 
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1.2 - AMBIENTE COMPUTACIONAL 


0 suporte computacional para o si sterna CARTA tem a se 
guinte configuracao: 



I9CK RK02 ( 2) RIC05(2) TV 10(2) CONTAL 


Fig. 1 - Configuracao do sistema CARTA. 

0 disco e de fabricacao da SYSTEMS INDUSTRIES, e consta de dois "drivers" 
para disco do tipo "Winchester", os quais emulam quatro discos RM02 da 
DEC, com capacidade total de 268 M "Bytes", com acesso dual (urn contro 
lador para PDP 11/34 e urn control ador para VAX lf/780). 

0 visual izador a ser utilizado (COHTAL) tera tela de 
512x512 pontos, com 8 bits para cada ponto, e sua capacidade de mos 
trar 256 tons de cinza devera ser bastante utllizada para a correta lo 
calizacao de pontos de controle. Este equipamento esta atualmente 1n£ 
talado em Cuiaba devendo ser transferido para Cachoeira Paulista para 
Integrar-se ao sistema. 

Tambem foi veriflcado (Mikhail and Paderes, 1983) que, 
para obter-se acuracia de menos de 1 "pixel" no processo, e impresci£ 
divel que os pontos de controle possam ser local izados com precisao 
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compatTvel. Uma futura expansao do sistema devera Incluir a compra de 
um processador matrlclal ("array processor"), o qual aumentarS cm 
to a rapldez de processamento do si sterna. 

2. DESCRICAO FUNCIONAL DO SISTE MA 

2.1 - DESCRICAO GLOBAL (ENTRADAS/SAlDAS) 

Do ponto de vista funcional, deve-se entender o CARTA co 
mo o sistema capaz de fornecrr, a partir de uma CCT e dos dados de efe 
merlde e de atitude do satelite, uma CCT da safai corrigido geometrica 
mente para uma projecao cartografica desejada. Se dispornveis, mapas 
com pontos de controle deverao ser Ntilizados; tambem podera ser gera 
da uma matriz de "break-points" (pontos de ancora) a ser utilizado no 
sistema VAX. De maneira esquematica, as entradas e saTdas serao: 

a) ENTRADAS 

1) CCT: Contem os dados de video em formato BIP2 (Secao 5.1). 

2) DADOS DO PCM: No caso do LANDSAT 1, 2 ou 3, tais dados coin 
tern contagens para "roll", "pitch" e "yaw", fornecidas a ca 
da segundo; para o LANDSAT 4, sao fomecidos os parametros 
de Euler a cada 4 segundos; as derivadas sao obtidas a cada 
0.5 segundo, a partir dos dados dos giroscopios e da posi^ 
cao e velocidade a cada 4 segundos. 

3) DADOS DA BFET: Fornecidos pelo programa Gylvax para os sa 
telites 1, 2 e 3, ou contidos nos dados de telemetria (Ban 
da S) para o satelite LANDSAT-4. 

4) FERRAMENTAS AUXILIARES: Utilizadas para a defirn'cao de GCP. 
Um exemplo seriam os mapas. 

5) ICT: Substitui os dados do PCM e das efemerides, e contem 
parametros para correcao, previamente gerados. A ICT e gera 
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da atualmente no DGI para o LANDSAT 1, 2 e 3 e sera um saT 
da opclonal deste si sterna para suportar o si sterna atual de 
geracao de fllme do DGI. 

b) SAIDAS 

1) CCT: Contem cena (ou mosaico) corrlglda geometricamente pa 

ra projegao especificada, com dados pertlnentes (Ex.: GCP 

utlllzados, nTvel de correcao, dados de efemerlde e PCM.ava 
llacao), no formato superestrutura. 

2) MATRIZ DE "BREAK- POINTS": A ser utlllzada para geracao de 
fllme corrigldo no EBIR, a partlr do dado original. 

3) ACI: Arqulvo com parametros de correcao de Imagem. 

4) "FLOPPY DISK": Contem um pedaco da Imagem corrlglda a ser 
utlllzado para entrada em sistemas que utillzem microcompu 
tadores. 

2.2 - FLUXOGRAMA DO SISTEMA CARTA 


0 fluxograma do CARTA consta de 8 modulos que executam 
as dlversas tarefas necessarias, a saber: 

a) MODULO 0 - GERENCIAMENTO: Envolve os procedlmentos de especlfi^ 
cacao das tarefas a ser reallzadas, bem como a verificacao dos 
procedlmentos ja efetuados e dos arqulvos resi dentes. Como o 
CARTA sera utlllzado em modo de producao, este modulo sera Iri 
teratlvo e de facll utlllzacao, permltlndo crTtica dos comaji 
dos do usuarlo; 

b) MODULO 1 - DETERMINACAO DOS PARAMETROS DA IMAGEM: Abrange todo 
o processo de equaclonamento da atltude e efemerldes do satelj^ 
te p Inclulndo calculo de "roll", "pitch" e "yaw". 
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c) MODULO 2 - AQUISICAO DE CCT: Abrange o processo de aqulslcao e 
geracao do arqulvo de Imagens a ser cor/lgldas. 

d) MODULO 3 - BIBLIOTECA DE GCP: Crlacao, gerenclamento, manuten 
cao e visual Izacao da blblloteca dos poitos de control e. Irt 
cltil tambem o fornecimento de GCPs da Imagem ao modulo de equ£ 
cion amen to. 

e) MODULO 4 - EQUACIONAMENTO GEOMETRICO DA IMAGEM: Determlnacao 

da matrlz de correcao da 1mag*m ("break-points") a partlr do 
conjunto de coordenadas e da velocldade do satellte, e dos da 
dos de "roll", "pitch", "yaw" e efemerldes, fornecldos pelo mo 
dulo 1. Caso exlstam GCPs, estes serao utlllzados para reflnar 
esta correcao. 

f) MODULO 5 - REAMOSTRAGEM: A partlr da matrlz de correcao obtem 
-se a Imagem LANDSAT corrlglda, devendo este modulo ser otlnri 
zado em termos de processamento. 

g) MODULO 6 - MOSAICO: Realiza a uniao de cenas LANDSAT para ob 
ter Imagem combined^ que contem area de Interesse ou obedece 
a grade padrao UTM. 

h) MODULO 7 - GERACAO CCT SAlDA: Fornece uma flta CCT em forma to 
especlflcado, com Imagem (ou mosalco) na projecao desejada. 

Para uma melhor visual Izacao, o fluxo de processamento 
do CARTA esta mostrado na Flgura 2 a seguir. 


I 


ITO 



Fig. 2 - Fluxo de processamento do 



si sterna CARTA. 





2.3 - DEFINICfiO DOS ARQUIVOS 


A def Inlcao dos arqulvos compart Hhados pelo si sterna con 
tempi a um conjunto de deflnlcoes a ser especlf Icadas a segulr. Inic1aJ_ 
mente sao esbocados os arqulvos de acordo com o esquema apresentado na 
Figure 3. 

Os mnemonlcos aprerantados na Flgura 3 denotam os seguln 

tes arqulvos: 

ICT - fita para correcao de Imagem (orb.), 

ACI - parametros para correcao (por cena), 

IMO - arqulvo com dados de Imagen (original), 

BPC - blblioteca de pontos de control e, 

BRP - arqulvo con "bre:\-po1nts" (pontos de ancora), 

PCI •• pontoc dc control e na Imagem (para cena), 

STA - Informacao sobre o slstema ("status file"), 

1HC - Imagem corrlglda, 

IMM - Imagem mosalco, 

DIRM - diretorio de mapas correspondentes e cena, 

AUX - dados auxll lares (anotacoes CCT), 

SAT - parametros de cada satellte. 
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Fig. 3 - Distribuicao das arquivos do slstema CARTA. 
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sentados: 


A seguir, serao descritos os arquivos anteriormente apre 


A.) BPC (Biblioteca de Pontos de Controle) 

Contem os pontos de controle separados em registros. Ca 
da GCP sera correspondente a uma local izacao em latitude e longitude, 
sendo guardado sempre o melhor GCP para uma dada coordenada. A compos^ 
cao de cada registro sera: 

- tTtulo do GCP, 

- identificacao do GCP, 

• latitude, longitude e altitude, 

- escala, folha e classe da carta correspondente, 

- origem e ordem do GCP (foto, mapa, 1? ou 2? ordem), 

- tipo do GCP: 1$ ordem (astr. , "doppler"), 2$ ordem, 

- data de insercao na biblioteca, 

- operador responsavel pel a insercao, 

- numero de utilizacoes, 

- merito acumulado (media, desvio padrao), 

- posicao dentro do "chip", 

- "chip" com o ponto (tamanho = 32 X 32), 

- banda do "chip". 
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B. ) PCI (Pontos de Control e na Imagem) 

Content os pontos de control e referentes a imagem a ser 
corrigida. Para cada registro ter-se-a: 

- identificacao do GCP, 

• latitude e longitude, 

- linha e "pixel" na imagem. 

C. ) STA (Arquivo de "status") 

Relaciona os servicos a ser executados, com todos os pro 
cessamentos realizados por servico, reservando urn registro para cada 
pedido, a saber: 

- ID do pedido. 

- ID das imagens correspondentes (content PCT). 

- Descricao servicos (produtos saTda, processamento): 

a) projecao desejada (UTM, SOM), 

b) geracao de filme no EBR, 

c) bandas a ser utilizadas, 

d) refinamento por GCP, 

e) tamanho de linha e "pixel" de saTda, 

f) escala do prodi'to final (caso de filme), 

g) informacoes para mosaico: 

- ID das imagens a ser mosaicadas. 
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- lat. , long, dos centros, 

- linha, “pixel" Inlclal para mosaico, 

1) tamanho da matrlz de “break-points". 

Processamentos reallzados (arquivos gerados) 
a) Bandas 

- faixas carregadas, 

- indice da imagem, 

- arquivos correspondentes: 

. entrada: 0 = existente, 

1 = nao-existente, 

. saTda : 0 = existente, 

1 = nao-existente, 

. mosaico: 0 = nao-existente, 

1 - exlstente, 

- correcoes realizadas: 

. "line-length": 0 = nao-inlciado, 

1 * iniciado, 

2 = terminado, 

. radiometrica: 0 = nao-iniciado, 

1 = Iniciado, 

2 * terminado. 


- reamostragem: 


. faixa da imagem de saTda: 0 = nao-iniciado, 

1 ■ iniciado, 

2 ■ terminado. 


- mosaico: 0 ■ nao-iniciado, 

1 « Iniciado, 

2 * terminado, 

. b) Arquivos auxiliares: 

. ACI: 0 = nao-criado, 

1 = criado, 

. PCI: 0 * nao-criado, 

1 * criado, 

. BRP: 0 = nao-criado 
1 = criado, 

c) CCTs realizadas. 

D. ) IMP (Arquivo de Imagem) 


Contem a imagem a ser processada no forma to BSQ (Bandas 
sequenciadas) tendo como entrada uma cena MSS. Para cada arquivo (uma 
banda) o conteudo sera: 

- “Header" de identificacao (urn bloco de 512 "bytes") com: ID 
cena, Banda, Faixas, TTtulo, Processamentos realizados e Ijn 
formacoes adicionais (ex.: "qualidade"). 
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- Dados da imagem com correcao do "line length". Ao final de 
cada linha serao colocados o "line-length" e o tempo corre± 
pondente ao "sweep". 

E. ) ACI (Arquivo para Correcao de Imagem) 

Gerado nos mesmos moldes da ICT, contera os parametros 
de entrada do "modelo rigido" para a cena, como a seguir: 

- "Header" com: 

. ID da cena, 

. PCT (tempo do centro da cena) designado por operador 
(ANO:DD:HH:MM:SS), 

. tempos siderais (ST1 - 0 HS., ST2-24 HS.), 

. numero de registros, 

. coeficientes do polinomio para atitude refinada, 

. tempo do primeiro grupo de dados, 

- M N" registros com: posicao (e suas derivadas), para cada se 
gundo (interpolados pelo modulo de determinacao de parame 
tros), alem dos dados de "roll", "pitch" e "yaw". 

F. ) BRP (Matriz de "break-points") 

Utilizada na reamostragem para geracao de imagem corrigi^ 

da com: 

- "Header" de identificapao que contem o ID da cena, indica 
dor de matriz direta/lnversa, numero de registros totais.nu 
mero de registros (NL) e numeros de "break- points"/ 1 inha 
(NBP). 
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- "NL" registros com LIN, PIX, LAT, LONG, e X.Y (na proj. 
UTM) com os "break-points" de cada linha. 

G. ) AUX (Dados Auxi Hares) 

Contem dados auxll lares, tais como: 

- anotacoes CCT de entrada (transferidos para a CCT de saida), 

- "tic marks" (references de locallzacao). 

H. ) SAT (Para metros de satelite) 

Contem os parametros de cada satelite utlllzados pelo mo 
delo de correcao. Oevem incluir: 

- Parametros relaclonados ao elipsoide: 

. semi-elxo maior (A), 

. semi-elxo menor (B), 

. achatamento (F), 

. orlgem no si sterna geocentrico (XO, YO, ZO), 

- parametros relaclonados ao espelho: 

. amplitude de movlmento (Rad.), 

. tempo util de varredura, 

. frequencla do movimento (Rad./seg.), 

. frequencia de varredur* (var./seg.), 

. constante de amortecimento, 

- parametros relaclonados com a Imagem: 

. numero de varredura/imagem, 
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. numero de llnhas/varredura, 

. numero de "pixel s"/l Inha), 

- parametros relacionados ao satelite: 

. altitude nominal , 

. resol ucao no sent i do da orbita, 

. resolucao no sentido transversal a orbita, 

. angulo de visa da por linha de varredura, 

- parametros relacionados com projecao. 

2.4 - POLlTICA DE ARQUIVOS 

A crlacao dos arquivos descritos anteriormente esta con 
dicionada por urn conjunto de diretrizes, onde algumas definicoes fi_ 
nais serao tomadas apenas durante o projeto detalhado. 

A) Arquivo de "status" (5TA) 

0 arquivo STA e a chave de acesso ao si sterna, pois rela 
ciona para cada pedido todos os arquivos associados. Como sera utiliza 
do urn arquivo para cada banda, no formato BSQ, o arquivo STA deve con 
ter a relacao dos arquivos associados, indicando sua utilizacao (origj_ 
nal, corrigida, mosaico). 

B) Arquivos de imagem (IMP, IMC, I MM) 

. Os arquivos de Imagem deverao ser pre-alocados na inicia 
lizacao do si sterna e, a cada arquivo, estara associado urn numero de 
Identificacao. Cada arquivo contera uma banda, no formato BSQ, e, no 
InTcio de cada tarefa, serao designados os arquivos para cada pedido. 
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Os arqulvos de Imagem terao parametros de tamanho flxo, 
Incluindo tamanho do "pixel", o numero de reglstros e o numero de li^ 
nhas. Deve-se notar que os arqulvos de imagem corrigida (KMC) terao ta 
manho diferente e serao descritos posteriormente. Estas dlferencas de 
vem levar em conta a mudanca de orientacao a ser imposta as imagens 
corrlgidas. 

3. DESCRICAO DO MODULO 0 - GERENCIAMENTO • 

, 3.1 - OBJETIVOS 

0 modulo de gerenciamento permitira ao operador do siste 
ma, a qualquer instante do processamento, uma visao completa sobre to 
dos os servicos em andamento. Para tanto, levara em conta as informa 
Coes constantes do arquivo de "status". 

A listagem do arquivo de "status" fornecera ao operador 
urn acompanhamento de todos os pedidos em andamento e dos procedimentos 
realizados em cada caso, indicando ainda os arquivos ja criados. 0 ope 
rador podera realizar ainda as fases de entrada de novos pedidos e ex 
clusao daqueles ja completados. 

3.2 - ENTRADAS/SAtPAS 

a) ENTRADAS: 

- formulario de pedido, fomecido pelo Servico de Atendimento 
ao Usuario do DGI; 

- arquivo de "status" (STA). 

b) SAIDAS: 

- "display" do STA (em impressora ou em terminal); 


- STA com pedido. 
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3.3 - AMB I ENTE/FLUXOGRAMA 


ORIGINAL PAL : . 

OF POOR QUALITY 



Fig. 4 - Ambiente para o modulo de Gerenciamento. 



-listagem 
-01 8 PLAY 


Fig. 5 - Fluxograma para o modulo de Gerenciamento. 
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4. OESCRICfiO DO MODULO 1 - DETERM I NAC AO DOS PARAMETROS DA I MAG EM 

4.1 -OBJETIVOS 

Este modulo fara a crlacao do arqulvo ACI a parti r dos 
dados de telemetria, recebidos da flta PCM e, no caso dos satelltes 
LANDSAT 1, 2 e 3, da flta de efemerldes (BFET). Se ja tlver sido gera 
da uma ICT para esta orblta, este modulo lera seus dados. 

4.2 - ENTRAP AS/SAlDAS 

a) ENTRADAS: 

- flta ICT da orblta, OU 

- flta PCM com dados de telemetria da orblta, E 

• flta de efemerldes, gerada pelo programa GYLVAX para LANDSAT 
1. 2, 3; 

- registro do STA (arquivo de "status") com descricao do servj 
CO a ser executado. 

b) SAIDAS: 

- arqulvo ACI com os parametros para correcao da imagem; 

- arqulvo com dados de telemetria, lidos da flta PCM, para uti_ 
llzacao no si sterna de geracao de Imagens MSS no VAX; 

- atuallzacao do registro correspondente no STA. 

4.3 - AM BIENTE/FLUXOGRAMA 

Duas observacoes sao pertlnentes: 1) a BFET nao e nece£ 
sarla para o LANDSAT 4, pols os dados de efemerldes estao incluTdos na 
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telemetrla (PCM). 2) As entradas podem ser substltuTdas por uma fi_ 
ta ICT. 



(PARA MSS DO VAX) 

Fig. 6 - Ambiente para o modulo de Determinacao dos Parametros. 

- DESCRICAO DOS DADOS DE ENTRADA 

Dependendo do satellte, os dados de trada serao dife 
rentes. No caso do LANDSAT 1, 2, 3, a entrada sera uma ICT, ou as 
tas PCM e BFET. Para o LANDSAT 4, a entrada e a fita PCM. 

A flta PCM para o LANDSAT 1, 2 e 3 content os errosdeatt 
tude para "roll", "pitch" e "yaw" a cada segundo, referenciados ao tem 
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po dado pelo reloglo do satellte, e a BFET contem as coordenadas carte 
slanas num slstema flxo com a Terra, a cada segundo, referencladas ao 
tempo GMT. A diferenca entre o tempo GMT e o reloglo do satellte £ da 
da por uma constante ("time-bias") que sera Informada pelo operador. A 
flta ICT contem os erros de atltude e a poslcao do satellte dados por 
latitude e longitude geodeslcas do ponto nadir e da altura do sateM 
te, referenclados ao tempo do satellte. 

Para o satellte LANDSAT 4 a flta PCM contem: 

1) A atltude do satellte em relacao a urn referenclal Inerclal com 
orlgem no centro da Terra, dada em quatemlos, a cada 4,096 
segundos. 

2) As lelturas dos glroscoplos (unldade de referenda Interior) a 
cada 0,512 segundos. Estes dados vao posslbllltar o calculo da 
atltude em relacao ao referenda! Inerclal dentro dos Intern 
los de 4,096 segundos. 

3) As coordenadas carteslanas do satellte em relacao ao referen 
. cla! Inerclal a cada 4,096 segundos. 

Para este ultimo satellte, a partlr das coordenadas do 
satellte e dos quatemlos, sao calculados os erros de atltude "roll", 
"pitch" e "yaw", a cada segundo, e as coordenadas do satellte no slste 
ma flxo com a Terra. Todos estes dados estao referenclados ao tempo do 
satellte. 

4.4 - DESCRICAO FUNCIONAL 


Inlcialmente, a funcao val lcr a Identlflcacao da Imagem 
do STA. Se for LANDSAT 1, 2 ou 3, pergunta-se se a entrada e ICT ou 
PCM ♦ BFET. Se for o segundo caso, pergunta-se o "time code" da ortrt 
ta. 


I 
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Fig. 7 - Fluxograma para o modulo de Det^rminacao 
dos parametros. 
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5. DESCRICAO DO MODULO 2 - AQUISICAO PE CCT 

5.1 - OBJETIVOS 

Este modulo e responsavel pels transference da imagom 
gravada em CCT para o disco do si sterna, que executa ao mesmo tempo cor 
recoes radlometricas e de "line-length". Os dados auxiliares contidos 
na CCT sao copiados para disco para posterior inclusao na CCT de saT 
da. 0 trabalho realizado por este modulo e sempre dirigido por urn ped^ 
do realizado por urn arquivo de pedidos. No caso de urn pedido de mosa2 
co, serao carregadas duas cenas. 

5.2 - ENTRADAS/SATDAS ‘ 

a) ENTRADAS: 

- CCT com imagem nao-corrigida (forma to BIP2), 

- registro do STA com espec if i cacao de servico. 

b) SAlDAS: 

- arquivo em disco no formato BSQ (urn arquivo/banda espectral), 

- arquivo com dados auxiliares (copiados da CCT). 

5.3 - AMBIENTE/FLUXOGRAMA 


a) AMBIENTE: 



\ 
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Fig. 8 - Anbiente para o modulo de Aquisigao. 
b) FLUXOGRAMA: 

A primeira fase deste modulo consiste na entrada da iden 
tificacao de urn pedido a ser satisfeito. No terminal do operador sao 

entao mostradas a identifi cacao da cena a ser carregada em disco, ou 

. 

das cenas no caso de mosaico, as bandas espectrais e as eventuais fa^ 
xas ja carregadas. Se todas as faixas ja foram carregadas, uma mensa 
gem e enviada ao operador e volta-se ao inTcio da funcao. 

• i 

. 

Apos a montagem da fita na unidade, urn teste e realizado 
para a confirmacao de que se trata da fita desejada. Em caso afirmatf 
vo, a fase de carga tern inTcio. A imagem e transferida para o disco 
faixa a faixa e o STA e atualizado indicando a presenga da faixa no ar 
quivo da imagem. Este procedimento e efetuado para as duas fitas que 
con tern uma cena e, em caso de contcccao de mosaico, duas cenas sao car 
regadas. Serao realizadas usual mente pre-processamentos de equal izacao 
radiometrica e corregao do "line length". 0 fluxograma apresentado na 
Figura 9 ilustra o processo. 




Fig. 9 
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• Fluxograma para o modulo de Aquisicao. 







6. DESCRICAO DO MODULO 3 - BIBLIOTECA DE GCP 


6.1 - OBJETIVOS 


Este modulo trata da construcao e manutengao de uma tn 
blloteca de pontos de controle (GCPs). A operagao desta biblioteca en^ 
volve as fases de entrada, recuperagao, exclusao e listagem de pontos 
de controle. A biblioteca e construTda de forma similar a disponTvel 
no Goddard Space Flight Center (Niblack, 1981). 

6.2 - ENTRADAS/SATDAS 

a) ENTRADAS: 

- arquivo de imagem (IMO), 

- cartas, 

- transparencias/imagem em papel , 

- medigoes geodesicas/astronomicas, 

- ACI. 

b) SAIDAS: 

- registros da biblioteca de GCPs (BPC): 

. operador, 

. local izagao geodesica, 

. escala da carta, 

. identifi cacao, 

. data de insercao, 

. folha, 

. Jtdos estatTsticos, 
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. numero de utilizacoes, 

. merito acumulado (qualidade), 

. posicao dentro do "chip", 

. "chip", 

. banda do chip, 

- mapa de distribuicao dos GCPs sobre uma dada area, 

- PCI (Pontos de control e da imagem de trabalho): 

. identificacao, 

. latitude e longitude, • • 

. linha e coluna, 

. altitude. 

6.3 - DESCRICAO FUNCIONAL 
A) IDENTIFICACAO DE GCPs 

Estes procedimentos compreendem todas as facilidades pa 
ra a identificacao dos pontos a ser considerados para insercao na bj^ 
blioteca, e inclui os seguintes passos: 

- selecao de GCPs utilizando cartas e transparencias, 

- caracterTsticas dos pontos: altura, latitude e longitude, 

- armazenamento das informacoes: altitude, atitude, identify 
cacao. 

' B) LOCALIZACAO DO PONTO NA IMAGEM 


Esta fase compreende as rotinas para relacionar o ponto 
na carta com o "pixel" na imagem, as quais fazem uso do equacionamento 
gecmetrico da imagem (Modulo 4). Os procedimentos incluem: 
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• estimativa da local izacao do ponto na imagem, 

- visual izacao da regiao da imagem que compreende o ponto com 
plena resol ucao, 

- local izacao refinada do ponto na Imagem: 

. "zoom" com reamostragem, 

. mudanca de banda espectral , 

. realce radiometrico, 

. principal componente, 

. local izacao manual do GCP na imagem, 

- armazenamento do "chip", 

- detecao de erros grosseiros na local izacao do GCP ("blunder 
detection"). 

• C) TESTE DE CONFIABILIDADE 


Estes procedimentos procuram medir a facilidade de corre 
ta identificacao de cada ponto a ser inserido na biblioteca. 0 ponto e 
comparado com seus vizinhos e sao efetuadas medidas estatisticas para 
determinar sua qualidade intrTnseca. 

D) RECUPERACAO DE GCP PARA CENA ESPECIFICADA 


Os procedimentos de recuperacao sao utilizados (durante 
o processamento de uma cena) para refinar o modelo de equacionamento 
geometrico. Tais rotinas envoi vem a recuperacao e identificacao dos 
GCPs contidos naquela cena, com a criacao do arquivo PCI (Pontos de 
Control e na Imagem), e incluem: 

- selecao dos GCPs apropriados em uma biblioteca, dada uma re 
giao. 


- estimativa de local izacao do ponto na imagem, i.e. (linha, 
"pixel") > f(atitude, orbita, prop, sensor., dados auxilia 
res), 

• local izacao refinada do ponto na imagem: modos manual e auto 
matico (SSDA e correlacao de bordas). 

- apresentacao de mapa de distribuicao dos GCPs na cena. 

E) GERACftO DE GCPs SUPLEMENTARES 

A geracao de GCPs suplementares envoi ve a local izacao de 
pontos diretamente a partir da imagem digital. Tal passo pode ser ne 
cessario em funcao do aumento de GCPs para uma dada cena. Este passo 
inclui as seguintes rotinas: 

- local izacao refinada do ponto na imagem (ver fase B), 

- util izacao dos parametros refinados de altitude para calc^ 
lar latitude e longitude, a partir da linha e "pixel", 

- armazenamento do "chip" e das informacoes pertinentes. 

F) EXCLUSAO DE GCPs NA BIBLIOTECA 

A exclusao de GCPs da biblioteca sera feita a partir de 
avaliacoes periodicas de qualidade e desempenho de cada ponto da bi^ 
blioteca. As funcoes necessarias sao: 

- selecao de GCPs apropriados na biblioteca, dada uma regiao, 

- apresentacao do mapa de distribuicao dos GCPs existentes na 
area de interesse, 

• exclusao manual dos pontos indesejaveis. 
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7. DESCRICAO DO MODULO 4 - EQUACIONAMENTO GEOMETRIC DA IMAGEM 

7.1 - OBJETIVOS 

Este modulo e responsavel pelo calculo do relacionamento 
entre a geometria estimada do satelite (altitude e efemerides) com a 
projecao cartografica e escala desejada. Este equacionamento e obtido 
a partir de urn modelo matematico que relaciona o vetor de visada do sa^ 
telite - relativo a urn par- linha, "pixel" - no sistema do satelite - 
com o correspondente num sistema inercial. Uma versao mais simplifies 
da deste modelo e encontrada em Forrest (1981) e uma abordagem comple 

ta em Mikhail e Paredes (1983). Neste caso,.a atitude do satelite e es^ 

tlmada por urn modelo parametrico, similar ao "modelo rTgido" descrito 
na proposta do Environmental Research Institute of Michigan (ERIM, 
1982), sendo a diferenqa mais significativa aquela do metodo de refina 
mento dos dados por meio de pontos de control e que pode fazer uso tan 

to de filtragem estocastica (Kalman), como de mTnimos quadrados. Meto 

dos mais sofisticados envolvem o uso de modelos com variaveis fisicas, 
nao-parametricas (Friedmann et alii, 1983). 

7.2 - ENTRADAS/SAlDAS 


As entradas do Modulo 4 serao os seguintes arquivos: A) 
ACI (arquivo com parametros de correcao da imagem); B) IMO (arquivo com 
dados de imagem); C) PCI (arquivo de pontos de controle referentes a 
imagem a ser corrigida). 

As saTdas serao os seguintes arquivos: A) BRP (matriz 
de pontos de ancora ("break-points") numa grade 9x9 que contem de£ 
vios delta X e delta Y; B) atualizacao do arquivo PCI que adiciona da 
dos referentes ao desempenho dos pontos de controle utilizados. C) "Tic 
marks" que indicarao as coordenadas geometricas da imagem corrigida; 
D) ACI, com a adicao dos coeficientes do polinomio que descreve a va 
riacao refinada de "roll", "pitch" e "yaw". 
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ORIGINAL PASS 13 
OF POOR QUALlVYj 

7.3 - AMBIENTE/FLUXOGRAMA 


AMBIENTE: 



Fig. 10 - Ambiente para o modulo de Equaclonamento. 

FLUXOGRAMA 

0 processamento deste modulo sera reallzado em duas fji 
ses. Caso estejam disponTveis para a imagem pontos de controle, os da 
dos de atitude ("roll", "pitch", "yaw") serao Inlcialmente refinados 
para aumentar a precisao do modelo. A fase final do processamento en 
volvera o calculo dos coeflclentes de distorcao para cada urn dos pon 
tos da matriz BRP. Serao descritas a seguir as fases do processo. 

A fase de reflnamento se faz necessaria devldo a baixa 
precisao dos sensores de atitude dos satelltes LANDSAT 1, 2 e 3. Tal 
precisao, de acordo com os dados fornet idos pel a NASA, e da ordem de 
1.0 mlllrradianos, e para obter acuracia da ordem de urn "pixel", estl 
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ma-se (Rlfman, 1973) que cada componente preclsa ser conhecldo com exa 
tldao de .01 mlllrradlanos. Deste modo, toma-se clara a importancia 
do processo de reflnamento para adlclonar Informacao no calculo final 
do model o de correcao. 

Uma caracterTstica importante de urn si sterna de correcao 
geometrlca e o procedlmento utlllzado para o reflnamento. Neste caso, 
deverao estar disponiveis para o usuario dols processos distintos: mT 
nlmos quadrados e flltragem estocastica. No primeiro caso, procura-se 
ajustar urn pollnomio ao conjunto de pontos de controle e calcula-se o 
resTduo obtido par.*» cada ponto. Numa segunda passada, sao omitidos os 
pontos onde foram observados desvios significativos. Este procedlmento 
e mals simples, embora funclone melhor com grande numero de pontos de 
controle. 

No caso de flltragem estocastica (Caron and Simon, 1975; 
Wong et alii, 1978), o refinamento faz uso de urn filtro de Kalman para 
atuallzar as estimatlvas iniciais dos coeficientes de atitude por meio 
de pontos de controle. Os coeficientes de atitude sao representados por 
urn vetor cujos componentes definem as realizacoes polinomlais das se 
rles temporals para "roll", "pitch" e "yaw". Mals expllcitamente, cada 
urn desses angulos tern sua variacao dentro da Imagem representada por 
urn pollnomio de tercelro grau, cujos coeficientes sao exatamente os 
elementos do vetor de estado do slstema. 

A estimacao dos coeficientes (vetor de estado) e feita 
de manelra sequencial: a estimativa a priori do vetor e progress ivamen 
te refinada tendo em vista as diferencas entre as local izacoes estima 
da e real para cada urn dos pontos de controle exlstentes. A formulacao 
do filtro pode ser resumida na equapao a seguir: 

X N0V0“ X ANT. + A(Y " Y REAL ^ * (1) 

Na equacao anterior a matriz de ganho A e determlnada a 
rartlr da formulacao do mTnlmo erro quadratlco medio (Kalman), e os 
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vetores Y representam as local Izacoes dos pontos de controle. Para um 
desempenho adequado do f11tro t a norma do erro residual para cada 6CP 
deve decrescer com o numero de observacoes reallzadas. 

0 calculo da local Izacao estlmada de cada 6CP e reallza 
do utlllzando os parametros do modelo obtldos pelo reflnamento ante 
rlor, a Figura 11 resume o processo. 



PARA CALCULO DE 
IIATMZ BRP 


Fig. 11 - Ilustracao do processo de reflnamento. 
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Apos o processo de reflnamento, ou caso nao exlstam GCPs, 
os desvlos para correcao da Imagem sao cal cul ados. 0 processo envoi ve 
o cal cul o, para todos os pontos da malha, da latitude e longitude de 
acordo com os dados de atltude e efemerldes do satillte. A segulr, pa 
ra cada ponto - de acordo com a projecao, escala e tamanho do "pixel" 
especlflcados - sao calculados os desvlos. 0 proredlmento do calculo 
de latitude e longitude envoi ve os segulntes passos: 

- determinacao da orlentagao do slstema do satelite em relacao ao 
slstema Inerclal ; 

- determinacao da reta de Interseccao com o ellpsoide (fornec^ 
coordenadas X, Y, Z no slstema Inerclal e usa os dados de "roll", 
"pitch" e "yaw"); 

- transformacao de coordenadas do slstema Inerclal para o slstema 
geocentrico de Greenwich; 

- determinacao de latitude e longitude para o ponto+; 

- determinacao das coordenadas X, Y na projecao desejada para o 
ponto. 

0 proximo passo envoi ve a determinacao dos desvlos a ser 
Impostos aos pontos de uma malha regularmente dlstrlbuTda em tomo da 
Imagem ("break-points"). Esta malha tera tamanho variavel, depenHendo 
da sua utlllzacao (geracao de fllme ou flta CCT) e do metodo utilizado 
para reamostragem (polinomlos ou Interpolacao linear por partes). 0 
calculo dos desvlos para cada ponto dependera da projecao, escala, ro 
tacao e tamanho do "pixel" desejadas, e esta resumldo na Flcura 12. 
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8. DESCRICAO DO MODULO 5 - RE AMO STR AG EM 

8.1 - OBJETIVO 

0 modulo de reamostragem reallza a parte custosa do pro 
cesso de correcao geometrlca. Aceltando como entrada uma matriz de pon 
tos de ancora (BRP), de tamanho varlavel , este modulo gera os parame 
tros necessarlos e reallza a correcao de cada ponto da Imagem de acor 
do com as especlflcacoes desejadas. 

8.2 - ENTRADAS/SAlDAS 

a) ENTRADAS: 

- arqulvo com Imagem a ser corrigida ("original"); 

- matriz de "break-points" (de tamanho varlavel); 

- arqulvo para Imagem de saTda ("corrigida"). 

b) SATDAS: 

• arqulvo com Imagem corrigida ("pixels" de 50 metros); 

• relatorlo com avallacao sobre Imagem. 

8.3 - AMBI ENTE/FLUXOGRAMA 

0 processo de reamostragem sera reallzado de manelra In 
versa, Isto e, dado urn ponto na Imagem de saTda procura-se local iza-lo 
na Imagem de entrada, Interpol ando a seguir para obter seu valor de 
clnza. Para determlnar os parametros do mapeamento inverso, a parti r. 
da matriz de pontos de ancora (saTda do modulo de equacionamento), sao 
propostos dols processos: 
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a) MAPEAMENTO POLINOMIAL: A partlr da matriz de. pontos da anco 
ra - que pode ser considerada uma lista de GCPs "sinteU 

cos" - e calculado um par de polinomios do 59 grau que fara 
o mapeamento para toda a imagem de saida. Os coeficientes dt_s 
te polinomic serao obtidos por um ajuste de mTnimos quadrados 
entre os pontos da malha, os qua is deverao obedecer a proje 
cao desejada. Este metodo e utilizado pelo ERIM e fazia parte 
da proposta apresentada (ERIM, 1982). Experimentos prelimina 


res realizados indicaram que possivelmente nem todos os 42 
coeficientes necessitam ser determinados, sendo suficiente 
equacoes do tipo: 

LINHA = a 0 + a x X + a 2 XY , (2) 

PIXEL = b 0 + b x X + .b 2 Y + b 3 XY + b^Y + b s XY 2 + 

b s Y 2 + b 7 Y 3 + b e XY 4 + b 9 Y 5 . (3) 

b) INTERPOLACAO LINEAR POR PARTES: Neste caso sera feito o calcu 
lo dos desvios para pontos de uma malha da imagem de saTda 
(correspondente a matriz de "break-points" invertida). Para 


cada ponto da malha de saTda, os desvios correspondentes se 
rao obtidos a partir de interpolacao linear em duas dimensoes, 
tendo em vista os pontos da malha de entrada mais proximos, 
que terao influencia inversamente proporcional a sua distan 
cia ao ponto considerado. 

Apos a determinacao da malha de saTda, os de c vios que 
indicam a posiqao do ponto desejado na imagem original serao cal cul a 
dos por interpolacao bilinear entre os quatro pontos da malha mais pro 
ximos, e sera entao necessario estimar o valor de seu nTvel d° cinza. 
(Moik, 1980). Isto pode ser feito lcncando mao de algoritmos de inter 
polacao. Neste sistema propoe-se o uso de duas alternativas: vizinho 
mais proximo (assume o valor de cinza no ponto menos distante) e con 
volucao cubica (o valor de cinza do ponto e a funcao de 16 vizinhos, 
conforme Bernstein (1976)). Este sistema pode ser resumido na Figura 13. 
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Fig. 13 - Ilustracao do processo de reamostragem 
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Deve-se notar que a imagem corrigida tera seus "pixels" 
com resolucao diferente da original (por exemplo, 50 metros) e orienta 
dos segundo direcao que pode nao coincidir com a imagem de entrada. No 
caso de mosaico, por exemplo, a cena corrigida estara sempre na orien 
tacao "norte-sul". Alem disso, metodos mais convenientes para o calcu 
1o da intensidade do ponto de saida poderao ser implementados future 
mente. Exemplos sao os metodos propostos: 

a) pelo ERIM para a compensacao das distorcoes inerentes ao seii 
sor MSS do satelite LANDSAT ("restoration") (Dye, 1975); 

b) pela McDonald, Detwiller and Associates (MDA) para a rotacao 
de imagens, levando em conta as condi goes de Nyquist (Friedmann, 
1981). 

9. DESCRICAO DO MODULO 6 - MOSAI CO 

9.1 - 0BJETIV0S 

0 modulo de mosaico esta destinado a fomecer, como pro 
duto, urn arquivo em disco que contem a combinacao de duas imagens ja 
corrigidas em produto padrao UTM (1,5 x 1,0 graus), ou pedido especial 
do usuario. 0 mosaico resulta na combinacao de duas cenas obtidas em 
datas distintas. Devido as facilidades de processamento do DGI - que 
incluem escolher o PCT para a CCT a ser gerada - e possTvel garantir 
a obtencao de urn mosaico correspondente a uma carta UTM apenas com duas 
cenas MSS do satelite LANDSAT. 

9.2 - ENTRADAS/SAlDAS 
• a) ENTRADAS 

- Arquivo de imagens que contem imagens corrigidas a ser contfn 
nadas. 
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b) SATDAi 

- A saTda deste modulo e urn arquivo de imagem, no mesmo forma 
to dos demais, com o resultado do processamento realizado. 

9.3 - AMBIENTE/FLUXOGRAMA 


Inicialmente, as duas imagens (ja corrigidas geometrica 
mente) serao submetidas a urn procedimento de equal izagao radiometrica 
que e considerado extremamente importante para garantir boa qualidade 
ao proauto final. A seguir, o operador devera deierminar no dispositi^ 
vo de visual izacao urn ponto de referencia comum as duas imagens. A se 
guir, por um processc de registro entre o mesmo ponto nas duas imagens, 
deti-.Tnina-se a regiao comum. 0 ponto de corte sera diferente para cada 
lima e sera escolhido de modo a minimizar o contraste radiometrico na 
regiao de juncao. A parti r da determinacao do ponto de corte para cada 
linha, sera realizada a operacao de "eu‘ ", justapondo os pedacos 
correspondentes das duas imagens no arquivo que contem o mosaico. Ad_^ 
cional mente, sera realizado um processo de filtragem passa-baixa 
(“smoothing") para diminuir a dife.renca entre os nTveis de cinza nos 
pontos de juncao (Mil gram, 1975). 

No tocante ao forma to de saTda, este devera consistir de 
uma imagem que corresponde a 1,5 graus de longitude por 1 grau de latj[ 
tude. Os "pixels" restantes serao preenchidos com zeros. E importante 
salientar que se consideram duas imagens como suficientes para obter o 
mosaico desejado, pois o INPE em sua estacao de geracao de CCT pode mo 
dificar o centro de qualquer cena. 

. 0 tamanho do "pixel" resultante devera ser dependente do 
formato a ser mantido na CCT/INPE. A principio, sugere-se que se consj^ 
derem duas opcoes: reamostrar o "pixel" para 50m, ou manter o tamanho 
original (57 x 79m). 

0 processamento do modulo de mosaico pode ser resumido 


• na Figura 14. 
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Fig. 14 - Ilustracao do modulo de mosaico. 

Desenvolvimentos futuros deverao incluir metodos mais so 
fisticados para o mosaico, que podera incluir procedimentos que modify 
quern o equacionamento da geometria da imagem, em funcao de urn melhor 
casamento na regiao de mosaico (Zobrist et alii, 1983). 

io. di :scricao do modulo i - geracao cct saIda 

10.1 - 0BJETIV0S 


- Este modulo e responsavel pela geracao de CCTs que am 
tern a cena geometricamente corrigida, armazenada em disco. 0 formato 
da CCT segue as normas da superestrutura recomendada pelo LTWG (ver do 
cumento INPE-2553/PRE/207). Este formato permite tanto bandas sequeii 



ciadas quanto bandas intercaladas por linhas (BSQ e BIL). Todas as in 
formacoes necessarias para a geracao da CCT devem constar do registro 
que con ten o pedido a ser satisfeito e demais arquivos do sistema. 



Fig. 15 - Ambiente para o modulo de geracao CCT saTda. 

10.2 - ENTRADAS/SAlDAS 
a) ENTRADAS: 

- arquivos com bandas corrigidas, 

- registro do pedido, 

- dados aux ilia res. 
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b) SAIDAS: 

- con junto de CCTs, 

- relatorio de execucao do modulo. 
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gloss Arid 


Este glossario content - para malor 1 egibil Idade -o si£ 
nlflcado de um conjunto de termos utlllzados ao longo do documento, os 
quals foram colocados em Ingles por se tratarem de indicacoes consagra 
das pelo uso. 

- BIP2: Formato de flta CCT em que sao atualmente distribuTdas 
as fltas MSS LANDSAT processadas pelo INPE. ("Band Interleaved 
by Pixel Pairs"). 

- BFET: Fita que content os dados estimados para efemerides do sa 
telite ("Best fit ephemerls tape"). 

- BSQ: Formato de fita compatTvel com computador onde cada banda 
da imagem do satellte esta armazenada separadamente das demais. 
("Band Sequential"). 

- BREAK-POINTS: Pontos de uma grade na imagem onde sao definidas 

as correcoes a ser executadas. 

- CCT: Fita compatTvel com computadur. 

- CHIP: Janela (pedaco) da imagem de tamanho fixo. 

- EPAS: Parametros de Euler que descrevem a atitude, fornecidos 
pelo satellte LANDSAT-4. 

- EBR: Gravador de fllme a felxe eletronico utilizado no DGI. 

- FLOPPY-DISK: Disco flexTvel e transportavel de pequena capaci^ 
dade, utilizado sobretudo em microcomputadores. 

- GCP: Pontos de referenda na terra cuja posicao geografica po 
de ser estimada ou conhecida. 
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GYRO: Glroscoplo inerclal que faz parte do slstema de mensur^ 
cao da atltude dlsponTvel no satellte LANDSAT-4. 

ICT: Flta com dados de correcao de Imagem, gerada atualmente 
no 


ID: Identlficacao da cena (unTvoca) DGI para correcao de lira 
gens LANDSAT. 

LINE-LENGTH: Comprimento da linha de imagem de satellte, que 
varla para cada varredura. 

PCM: Forma to em que esta a flta analogica que contem os dados 
de atltude do satellte ("pulse coded modulation"). 

PCT: Tempo do ponto central da cena LANDSAT ("picture center 
time"). 

PIXEL: Elemento mTnlmo constituinte de uma Imagem digital. 

ROLL , PITCH , YAW: Movlmentos nos tres eixos que compoem o sl£ 
tema de referenda do satellte (rolamento, arfagem e gulnada). 

SWEEP: Varredura correspondente a cada linha da imagem LANDSAT. 

TIC-MARKS: References para local Iza cao da Imagem corrlgida, 
que corresponderii a Indlcacoes de latitude e longitude. 

UJM: Projecao cartograflca (Universa. transverse Mercator) na 
qual sao construTdas cartas mllltaros. 


